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As it was shown recently, parametric resonance during preheating opens new possibili- 
ties for electroweak baryogenesis [1, 2] which can take place both during [3] and after [4] the 
resonance. While baryogenesis mechanisms described in [3] and [4] can appear in the same 
theory and even overlap in the course of the same parametric resonance, they are based on 
different dynamic effects. In the present paper we compare baryoproducing efficiency of both 
mechanisms and discuss underlying dynamical processes. The key difference between these 
mechanisms is that they are related to different time scales existing in the model. Baryogen- 
esis at preheating exists on rather short parametric resonance time scale, while baryogenesis 
after preheating is related to inflaton thermalization time scale, that can become very large, 
at least in classical dynamics [6]. The third time scale is the Higgs field thermalization 
time which generally lies between the two previous time scales. This intermediate time scale 
defines the freeze-out time of post-resonant sphaleron transitions and is important because 
sphaleron transitions may wash out the generated asymmetry if they'll keep going after the 
end of baryoproduction. While the wash-out processes as well as other fermion backreaction 
effects are beyond the scope of our numerical simulations, the very existence of intensive 
sphaleron transitions on the time scale considerably larger than baryoproduction time is a 
signal for possible problems. 

The resonance- localized baryogenesis [3] can rapidly produce considerable shift in topo- 
logical number. Unlike baryogenesis at first-order phase transition, it occurs during extended 
period of time on the resonance time scale, thus making the wash-out problem a bit easier 
to solve. However, in this case wash-out suppression requires rather short Higgs thermal- 
ization time which should be close to resonance time scale. An attempt to achieve this in 
(l+l)-dimensional theory being simulated in [3] leads to unacceptably high reheating tem- 
peratures T rh ~ 0.3 v where sphaleron transitions do not freeze out at all, see Fig. 8 of [3]. 
Even though this conjecture isn't directly applicable to realistic (3+l)-dimensional case, it 
shows that more efficient ways to prevent wash-out are highly desirable. 

A promising approach was suggested in paper [4] where it was shown that for certain 
parameter values baryogenesis keeps going much longer than the resonance itself. In fact, CP 
violation that drives baryoproduction is present here for very long time, and baryoproduction 
keeps going as long as the sphaleron transitions occur. This scenario is based on difference 
between inflaton and Higgs thermalization times. In our simulations Higgs field thermalizes 
much faster than inflaton field because it couples to gauge field and has nonlinear self- 
potential, which results in intensive rescattering of Higgs spectral modes. As we will show 
below, the presence of non-thermalized oscillating inflaton field can generate effective CP 
violation, which exists on the inflaton thermalization time scale, the longest one in the 
model. Of course, this CP violation can change the baryon number only via sphaleron 
transitions, so the baryogenesis actually occurs on Higgs thermalization time scale, freezing 
out together with the sphaleron transitions. The major advantage of approach [4] is that 
the baryoproduction-driving effective chemical potential exists here long after the sphaleron 
transitions vanish. In other words, it gives us a toy model of electroweak baryogenesis which 
is completely free from wash-out. 

An interesting peculiarity of baryogenesis at and after preheating is that both these 
scenarios can realize simultaneously, see Fig. 2. Once for the same initial parameters they 
generate baryonic asymmetry of opposite signs, reflecting the fact that underlying dynamical 
mechanisms are entirely different, the initial conditions determine not only the magnitude, 
but also the sign of the net shift in baryonic number. 
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1 Baryogenesis at preheating 



Dynamical principles of baryogenesis at preheating are relatively simple (for detailed de- 
scription of physical framework, Lagrangian and dimensionless units used below see Ref. [3]). 
Energy transfer from inflaton field into low-momentum modes of Higgs field results in high 
rate of inequilibrium sphaleron transitions that appear very soon after the beginning of res- 
onance (see Figs. 8 and 14 of paper [3]). On the other hand, in the course of resonance the 
Higgs v.e.v. (0*0) comes from zero to its post-resonant value close to 1. Similar to earlier 
models of baryogenesis at first-order phase transition (see e.g. [5]), the time evolution of 
Higgs field gives rise to effective chemical potential /i e s that appears through CP-breaking 
term — k|0| 2 e^F^ in Lagrangian. In the presence of sphaleron transitions, this chemical 
potential drives the Chern-Simons number N cs — ^ / A x dx in certain direction (see Figs. 11 
and 16 of Ref. [3] and Fig. 1 below), thus producing the fermions. In our case, the chemical 
potential is [3, 4] 

/Xeff = -4717t<9 O (0*0) (1) 

and has nonzero value only while (0*0) keeps changing. 

The transition of (0*0) to its stationary value isn't instant (see Fig. 10 of [3]); it takes 
about the same time as the resonance itself, which is generally much slower than the first- 
order phase transition. This is an important advantage of baryogenesis at preheating over 
previous scenarios, because extended baryoproduction time is essential for reducing possible 
destruction of the created baryons by post-resonant sphaleron transitions, although sphaleron 
freeze-out is related to Higgs thermalization and occurs on different time scale. 

A common way of estimating the baryoproduction in this and similar models is the use 
of diffusion equation for Chern-Simons number: 

^CS _ _ r /feff / x 

dt ~ sph T cff [Z) 

where r sp h is the rate of sphaleron transitions and T e g — certain effective temperature. Note 
that Eq. 16 of [3] coincides with Eq. (2) when T B = (the latter condition reflects the lack of 
fermionic backreaction in our simulations). Once the (0*0) value is increasing in the course 
of resonance, Eqs. (1) and (2) give the sign of AiV cs : 

oc oc <9 O (0*0) > (3) 

K K 

Eq. (2) can be used for rough estimate of total generated asymmetry by assuming final 
(0*0) to be 1: 

Aiv cs ~ ^(r sph ) (4) 

giving AiV cs ~ —20 for Fig. 1 and A7V CS ~ —30 for Fig. 2. While the first estimate is correct 
within factor of 2, the latter one has the wrong sign and is off by an order of magnitude. 
This descripancy cannot be explained by simplifications used for getting Eqs. (3) and (4). 

Indeed, it isn't obvious that the diffusion equation Eq. (2) can be reliably applicated to 
stronly nonequilibrium dynamics of parametric resonance. We check this point directly by 
substituting (0*(t)0(t)) and Y{t) measured in numerical simulations [3] into Eqs. (1) and 
(2), taking T e g = T r h and numerically integrating over time, see Figs. 1 and 2. The only 
source of uncertainty in this check is the effective temperature T e g present in Eq. (2). Exact 



3 



physical sense of this quantity isn't completely clear; however, Figs. 6 and 13 of Ref. [3] show 
that various effective temperatures stay rather close to final equilibrium value T rh and their 
time variations wouldn't considerably affect the time integral of Eq. (2). 

As is clear from Fig. 1, even for high reheating temperatures the Eq. (2) is adequate only 
for rough qualitative analysis of N cs evolution. It cannot reproduce correclty the short-time 
dynamics of N cs and provides no explanation for the end of baryoproduction at time ~ 1800. 
Both problems are unrelated to possible variations in T cS which are slow and disappear at 
time ~ 1500, see Fig. 6 of Ref. [3]. For lower T rh ~ 0.09-E sp h the diffusion equation (2) is 
no longer able to give correct sign of final asymmetry, see Fig. 2. It is natural to suppose 
that the evolution of N cs cannot be reduced to simple diffusion described by Eq. (2) and is 
sensitive to more subtle dynamical effects. A good example of complicated dynamics beyond 
the scope of Eq. (2) is the baryogenesis after preheating. 

2 Baryogenesis after preheating 

The driving force for baryogenesis after preheating [4] is the inflaton field that keeps co- 
herently oscillating after incomplete parametric resonance (see Figs. 2 and 3 of Ref. [4]). 
Inflaton zero mode decays through thermalization which in our case is very slow. Its oscil- 
lations considerably modify Higgs field effective potential 

V(0) = ^(|0| 2 -, 2 ) 2 + ^V|0| 2 (5) 

thus making inflaton zero mode (a 2 ) and Higgs v.e.v. (0*0) to oscillate coherently. In 
its turn, (0*0) oscillations generate periodic CP-nonconservation through Eq. (1) and also 
trigger the sphaleron transitions that become synchronized to inflaton/Higgs oscillations. 
Interplay between these two strongly correlated processes results into steady shift of N cs in 
direction opposite to that one predicted by diffusion equation (2). 

In this section we will analyze dynamics of Chern-Simons number N cs in the presence 
of oscillating ji c ^. Then, taking into account time correlations to Higgs-triggered sphaleron 
transitions, we show why iV cs actually shifts in direction opposite to baryogenesis at pre- 
heating. 

2.1 CP-(0*0) correlations 

Due to time derivative in /i e g, Eq. (1), periodical oscillations of (0*0) clearly affect N cs 
behavior. For k — 0, 7V CS is oscillating uncorrelated to (0*0) with its own characteristic 
frequency ujq = My/ ~ 1 (these oscillations are averaged out on the plots), while introduction 
of non-zero k makes the N cs dynamics dominated by n e $ even for relatively small k ~ 0.1, 
see Fig. 3. The amplitude of N cs oscillations at non-zero k can be estimated using equation 
of motion [3, 4] for A 1 field. Let's assume that N cs is oscillating around its steady integer 
value Nq which is the topological index of current system vacuum. Making large gauge 
transformation to vacuum with zero topological index and denoting the new Chern-Simons 
number as SN CS = N cs — N , one gets 

5N CS + 2(0*0)<5iV C s = • 2/^(0*0) (6) 
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(we assume the Higgs field to be homogeneous, i.e. neglect all its higher-momentum modes 
except the zero mode), and 



SN CS = -^— 2 ■ ^5(0*0) (7) 

Here £(0*0) = Ae tuJat (A is the amplitude of (0*0) oscillations), uj 2 = 4g 2 (0*0) is the 
frequency of inflaton oscillations, and ujq 2 = 2(0*0) = My/ 2 is the proper frequency of free 
N cs oscillations not visible on Fig. 3. Once g 2 <C 1, we can neglect uj 2 <C Uq 2 and further 
simplify Eq. (6): 

= /tL <9 O (0*0) L /X eff 

CS 2tt (0*0) 4tt 2 (0*0) {) 

Both eqs. (7) and (8) result in the phase shift between 8N C s and 5(0*0) being exactly 
| sign k. This phase relation actually determines the sign of baryonic asymmetry generated 
by T — CP correlations. 

2.2 r sp h-(0*0) correlations and baryoproduction 

Once (0*0) is oscillating periodically, time averaging of fi c g (1) and SN CS (8) will give zero. 
In our case this does not mean zero baryoproduction, because the sphaleron transitions occur 
synchronously to the same oscillations of (0*0), thus favouring certain values of /z e fr- 

The correlations between sphaleron transitions and (0*0) oscillations are discussed in 
more detail in paper [4]. It is quite natural, however, that most transitions occur at (0*0) 
close to its minimum (see Figs. 4 and 6 of Ref. [4]), once this corresponds to minimal E sp ^. 
Note that at these moments /i e fj oc <9 O (0*0) = 0, so (/x e ffF sph ) = 0. In other words, Eq. (2) 
fails even after accounting for correlations between /z e fr and r sph . Therefore, baryogenesis 
after preheating involves dynamical effects that are time-nonlocal. Once most topological 
transitions occur when (0*0) is at minimum, (0*0) value is decreasing during the whole 
half-period of (0*0) oscillation immediately preceeding the transition, so (9o(0*0) < and 
yUefr has non-zero value of certain sign. According to Eq. (8), 5Nqs/k oc <9o(0*0) < over 
the same half-period, see Fig. 3. Therefore, the Chern-Simons number is always shifted 
in certain direction for some period of time just before the sphaleron transitions, making 
transitions with AN CS = — sign k more probable than transitions in opposite direction. 

This description provides a simple estimate for net asymmetry created by baryogenesis 
after preheating. Assuming that every topological transitions effectively freezes the maximal 
positive or negative SN CS that existed 1/4-period before, one obtains from Eq. (7) (here we 
also use Eqs. 2.1 and 2.2 of Ref. 



dN cs ■ ( XAT \ r KgvL /a max y 

- sl s n K ■ ( 5iV cs) max r sph = — — J r sph (9) 



dt 

and from Eq. (8): 



ANcs -^ O (0*0)<O. (10) 



K 

(compare to Eq. (3)). Substituting into Eq. (9) parameter values from [3, 4], one obtains 



dN, 



dt 



' s ~-0.3«T sph (11) 
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in reasonable agreement with numerical results, see Fig. 4. 

Although baryogenesis after preheating is suppressed by inflaton-Higgs coupling g, it 
keeps going as long as topological transitions occur, see Eqs. (9) and (11), and isn't affected 
by the wash-out problem. It also shows how strongly the nontrivial dynamics can affect 
simple estimates in the style of Eq. (2). 

3 Conclusions 

Here we have shown that dynamical effects do play an important role in baryogenesis at / after 
preheating, what is typical for complicated nonequilibrium processes [7, 8]. A simple (1+1)- 
dimensional model of electroweak preheating-related baryogenesis, depending on initial pa- 
rameter values, may subsequently pass through two dynamically distinct stages with differ- 
ent baryogenesis mechanisms producing baryonic asymmetry of opposite signs. The second 
stage that occurs after preheating is completely beyond the reach of straightforward diffusion 
approach commonly used for analysis of electroweak baryogenesis scenarios. 
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Figure 1: (upper plot) Time evolution of Chern-Simons number for k = — 1 and high T r h = 0.33 v (solid line, 
identical to Fig. 11 of Ref. [3]) compared to prediction of Eq. (2) (dotted line). 

Figure 2: (lower plot) Same for reduced-energy runs of [3] with k = — 1 and T r h = 0.094 v. On this 
plot {Ncs(t)), solid line, differs from Fig. 16 of Ref. [3] due to larger statistics. Strong descripancy between 
observed (Ncs(t)} and that one predicted by diffusion equation (2) means that the baryogenesis at preheating, 
roughly described by Eq. (2), is followed here by some other baryoproduction mechanism — baryogenesis 

after preheating that moves N cs in opposite direction. 
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Figure 3: Comparison of phase relation between N cs and ((/>*(/>) oscillations for n = 0.1 (upper plot) and 
k = —0.25 (lower plot) shows that their relative phase shift is sensitive only to sign k. Both plots also confirm 
Eq. (8). The sphaleron transitions occur when {<j)*(f>} is at minimum and N cs is close to its mean value -/V , 
see text. Note that curves presented here are averaged over large ensemble of independent runs, so -/V has 

non-integer values —0.8 and 2.9, respectively. 
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Figure 4: Eq. (11) adequately describes dynamics of Chern-Simons number both for n = 0.1 (upper plot) 
and k = —0.25 (lower plot), demonstrating the lack of wash-out for baryogenesis after preheating. Unlike 
Fig. 2, both (Nqs) curves can be smoothly interpolated to origin, due to the suppression of baryogenesis 

during preheating at very low T r h ~ 10~ 2 v. 
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